Exact Scaling Functions for Self- Avoiding Loops and Branched Polymers 
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It is shown that a recently conjectured form for the critical scaling function for planar self- 
avoiding polygons weighted by their perimeter and area also follows from an exact renormalization 
group flow into the branched polymer problem, combined with the dimensional reduction arguments 
of Parisi and Sourlas. The result is generalized to higher-order multicritical points, yielding exact 
values for all their critical exponents and exact forms for the associated scaling functions. 
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In two dimensions, exact results for critical exponents, 
which describe power law dependence on a single relevant 
variable close to a critical point, are commonplace. By 
contrast, very few examples are known of exact scaling 
functions which depend on combinations of more than 
one such variable. In recent years there has been consid- 
erable progress in problems involving counting various 
restricted classes of random self-avoiding polygons pi. 
However, these are essentially one-dimensional in nature, 
and so far no rigorous results exist for the unrestricted 
case. Theoretically this latter case is perhaps more in- 
teresting because, in the scaling limit, it corresponds to 
an isotropic field theory, the n — > limit of the O(n) 
model. While much exact information is known about 
such critical theories in two dimensions, up to now no 
exact, nontrivial, scaling functions of more than one in- 
tensive thermodynamic variable (such as the equation of 
state) have been found. 

Recently Richard, Guttmann and Jensen 0] (here- 
inafter referred to as RGJ) have conjectured the ex- 
act form of such a scaling function for unrestricted self- 
avoiding polygons in the plane. In the ensemble in which 
each link of the polygon is weighted with fugacity x, this 
problem exhibits a critical point at some value x = x c . 
As x — ► x c from below, the mean perimeter (N), mean 
square radius of gyration (R 2 ) and mean area (A) all di- 
verge. Such self-avoiding loops provide a simple model 
for two-dimensional vesicles H, and in that context it 
is natural to weight the ensemble according to the area 
of each loop, thus defining the generating function for 
rooted loops 

GM(,, 9 ) = ^ P i;|/e-^ 



N,A 
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where p N \ is the number of such loops of a given perime- 
ter and area which pass through a given link of the lattice, 
and g (= — laq in the notation of RGJ) is the pressure 
difference across the vesicle wall, in units of kT. From 
the point of view of critical phenomena, this ensemble 
has a multicritical point at X = x c , g = 0, and in its 
neighborhood one expects H that the singular part of G 
has the scaling form (in the notation of RGJ) 



G%{x,g) = g e F({x c -x)g- 



(1) 



where 9 and </> are related to the conventional exponents 
v and a by 9/4> = 1 — a and <p = 1/2^, based on the as- 
sumptions that, at g = 0, the singular part of the mean 
number of such rooted loops behaves like (x c — a;) 1 ^" and 
the mean area (A) ~ {R 2 } ~ (x c -x)~ 2u ~ {N) 2u . In ad- 
dition a is related to v by hyperscaling: a — 2—2v. There 
is ample evidence from enumeration and other methods 
[ 1 1 to support these assumptions and the theoretical value 
[5| v = |, so that 9 = | and <j> = |. 

RGJ, in analogy with similar but analytically tractable 
enumeration problems |l||, assume that, as function of a; 
and q = c~ s , G^ satisfies some g-algebraic functional 
equation of finite degree. Together with the assumed 
values for 9 and </>, this leads, in the limit q — ► 1, to 
a Riccati equation for the scaling function F(s), whose 
solution is 



F(s)=b — lnAi(M) 
as 



(2) 



where Ai(#) oc J e ixt+it /3^ j g ^g Airy function, 
and 6o, b\ are non- universal constants. (g) deter- 
mines exactly, for example, the universal moment ratios 
(A P )/{A) P asi-> x c , and RGJ produce convincing evi- 
dence, based on extensive enumerations, that these pre- 
dictions are indeed correct. 

In this Letter, it is pointed out that (0) also fol- 
lows from a completely different argument, which invokes 
the physical reasoning of Ref. fl to relate this prob- 
lem to that of branched polymers, combined with the 
dimensional reduction arguments of Parisi and Sourlas 
H which map this latter problem to that of the Yang- 
Lee edge singularity in two fewer dimensions. From this 
point of view, the Airy integral then arises as the scaling 
limit of the partition function of the Yang-Lee problem 
in zero dimensions. Moreover, in this approach, the val- 
ues of the exponents 9 and <f> emerge without any further 
assumptions. 

From this perspective it is simple to generalize the 
conjecture of RGJ to higher-order multicritical points of 



self- avoiding loops with k relevant renormalization group 
(RG) scaling variables Vj . These may presumably be re- 
alized by tuning to critical values many-body interactions 
between nearby portions of the loop. When this ensem- 
ble is, in addition, weighted by the area of the loops, the 
generalization of (nT) to arbitrary k is 



G^ = 9 Bk F k {v x g 



-1/1/2 
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where yj is the RG eigenvalue of Vj . It will be argued that 
the exact values for these, at the fcth order multicritical 
point, are 



y s (k)=2{k-j + 2)/(k + 2) 



(4) 



and that the exact form for F^ is given in terms of a 
generalized Airy integral ^ c e~ v ^^ 9 dp where V(tp) — 

52 7=1 v i^ ~ 4> k+2 - The values given in (|j) agree with 
those derived from a generalized Flory argument, applied 
to an ensemble in which the first k renormalized virial 
coefficients vanish. That they should be exact in two 
dimensions was suggested earlier by Saleur M on the ba- 
sis of a postulated N = 2 supersymmetry. Here it is 
seen that they follow from the mapping to a simple zero- 
dimensional problem. However, our results, like those of 
Saleur, display a paradox in that k — 2, the obvious can- 
didate for the 6-point, yields values for the exponents 
which disagree with those of an exactly solvable model 
|19fl and with extensive numerical results. This is dis- 
cussed in detail later. 

Finally, in the generalized Airy integral it is possi- 
ble to take the limit g — * 0, thus recovering the scaling 
function in the original ensemble with no area-weighting. 
This comes from the appropriate saddle-point of V(ip), 
and therefore amounts to finding the root of a polyno- 

(r) 

mial. For example, for k = 2 it is found that G^-' = 

cqv 2 ' $(cit>i/v 2 ' ), where Co and c\ are non-universal 
constants, and the exact scaling function, for v<i > 0, is 

$>(*) = (a + (s 2 - I) 1 ' 2 ) 1 ' 3 + (s- (s 2 - l) 1 / 2 )^ 3 (5) 

where the branch cuts of the fractional powers are taken 
to lie along the negative real axis. 

We now give more details of the reasoning leading to 
these results, first discussing the case k = 1 considered 
by RGJ. The physical part of the argument is to regard 
the model with g > as presenting a crossover phe- 
nomenon: the negative pressure causes the vesicles to 
try to minimize their area, but there is competition be- 
tween this and the need to maximize the perimeter as 
x — ■> x c . Clearly for large g at fixed x < x c the vesi- 
cles should collapse into double- walled, branched struc- 
tures, but assume, as suggested by the numerical work 
of Ref. H that, at large enough distance scales, this will 
also happen as x — > x c for any fixed g > 0, consistent 
with the idea that there is an RG flow from the fixed 



point describing self-avoiding loops to that correspond- 
ing to branched structures. Thus (Q) has the form of a 
crossover scaling function M. We also assume that the 
structures which result are in the same universality class 
as conventional branched polymers (lattice animals with 
no cycles), in which all trees with the same total length 
are weighted equally. The theory of crossover scaling 
then asserts that the scaling function F(s) in (pi) should 
have a singularity of the form (s — s*) 1 ~ aBP , where asp is 
the entropic exponent for branched polymers. From this 
hypothesis various interesting results follow, for example 
that as g — ► the branched polymer singularity should 
occur at x = x c (g) = x c + s*g^ + ■ ■ •, which has been con- 
firmed in enumeration studies [g|, as well as various pre- 
dictions for the g-dependence of the critical amplitudes. 
In general, however, the functional form of a crossover 
function is very difficult to calculate, since the scaling 
variables at the new fixed point bear a complicated re- 
lationship to the original ones, which requires following 
the RG flow in detail. However in this example there are 
considerable simplifications. 

First state the problem in field-theoretic language, by 
writing the area of a given loop as 



A = 



G\ a {ri - r 2 )Jx(ri)J a (r 2 )d 2 r 1 d 2 r2 



(6) 



where J\ is the density of a current of unit strength flow- 
ing around the loop, and G\ a is the Green function for 
a U(l) gauge field A. (g) expresses the well-known fact 
that, in a two-dimensional gauge theory, the expectation 
value of a Wilson loop obeys a strict area law. In Ref. |J] 
it was used to compute the mean area of self-avoiding 
loops. The generating function G^ for rooted loops is 
the derivative with respect to the fugacity x of 



Z = <C- 9A >SAL = 



-y/5 f J\A*(Pr\ 

/SAL,A 



where the average is taken over self- avoiding loops (SAL), 
each weighted by x^, and over the gauge field, with the 
usual weight exp(— j j F Xa F\ a d 2 r). Self-avoiding loops 
may be mapped, in the standard way, to the n — > limit 
of an O(n) model. In this case it is useful to consider 
complex O(n) lattice spins s(r), so that the U(l) current 
J\ is the lattice version of (1/2i)(s*-c?as — cc), and the 
weights are Ilnn(l + x(s* (r)-s(r') + cc.)). 

The first observation is that, at n = 0, there are no vac- 
uum corrections to the gauge field propagator G\ a (as in 
the 'quenched' approximation in lattice gauge theories), 
so that, since A couples to a conserved current, the gauge 
coupling g is not renormalized. Its RG equation is sim- 
ply dg/d£ = 2g to all orders, so that it flows to infinity, 
where the irrelevant variable g~ 2 has RG eigenvalue —2. 
The other simplification is that, in the limit where g is 
large, the total length of the branched polymer is one 
half that of the perimeter of the loop, apart from correc- 
tions o(N). Thus the fugacity variable for the branched 



polymer problem, close to x c , is linearly related to the 
original fugacity x. 

Now recall the formulation of the branched polymer 
problem in d dimensions, given by Parisi and Sourlas 
H . This is the n — ► limit |10| of a theory of fields ip a 
(a = 1, . . . , n), weighted by e^^ where 



J a p>l ab 



cfir (7) 



where u p is the fugacity for p branches to meet at a given 
point, and v > represents self-avoidance. After shifting 
the fields to eliminate the ip 2 term, rescaling, and retain- 
ing only the most relevant couplings, the action takes the 
form 



S = 



I (E(5(^a) 2 + Vtya)) + AY,Mb)d d r (8) 



where V(ip) = tip — hip 3 + 0(ip 4 ). This theory is crit- 
ical at some value t — > t c +. Parisi and Sourlas |6[ ar- 
gued that, at n = 0, (jg) is equivalent to a supersym- 
metric theory. We follow the more direct transformation 
of Ref . 1 1 1 1 : define new combinations of the fields ip = 

i(Vl + (n-l)- 1 E2^)^^A- 1 (V'l-("-l)- 1 E2V'a), 

together with n — 2 other fields Xa (a = 3, . . . , n) which 
are linear combinations of (ip2, ■ • • ,ipn) orthogonal to 
^2^. Discarding terms higher than quadratic order 
in u> and the Xa (which may be shown to be irrelevant), 
the action has the form, at n = 0. 



S=±J(uj(-V 2 ^ + V'(i/j))-uj 2 

+£xa(-v 2 + v ,, WO) X a)<ri- (9) 



The integral over the n — 2 commuting fields Xa yields 
det(— V 2 + y")~(™~ 2 )/ 2 and so they may be replaced at 
n = by two anticommuting fields x an d X- The super- 
symmetry is made explicit by introducing anticommuting 
coordinates (9, 9) and a superfield \I> = ip + \ {Ox + &X) — 
j89lu, whence S may be written 



5' 



1 
A 



*(- 



-vis)*- 



V(*)) dfirdBdB (10) 



where V| 5 = V 2 +4d 2 /d9d9. This exhibits supersymme- 
try under rotations which leave r 2 + 99 invariant. Parisi 
and Sourlas o argued that this theory exhibits a remark- 
able property of dimensionality reduction (for a nonper- 
turbative proof see Ref. pl|): correlation functions whose 
arguments are restricted to ad- 2-dimensional subspace 
are the same as those for a non-supersymmetric theory 
in d — 2 dimensions, whose action is 



5- 



YL 



1 

A 



(iV(-v 2 )v^ + v(V0)) d d -\ (11) 



where in this case V(ip) = tip — \ip z ■ There is one sub- 
tlety: before dropping the irrelevant terms, the contour in 
-0 should be rotated, in this case parallel to the imaginary 
axis, so as to make the integral defined nonperturbatively. 
The potential therefore becomes ittp + ^iip 3 - Thus dO) 
is just the action for the scaling theory of the Yang-Lee 
edge singularity |Q, as discussed by Fisher [0|. From 
( |l0| ) it is seen that A has dimension (length) - and this 
is not affected by loop corrections, otherwise supersym- 
metry would be broken Jl5[ . It flows to infinity under the 
RG, and A -1 is irrelevant. However, it is a classic exam- 
ple of a dangerously irrelevant variable: it cannot be set 
equal to zero in the scaling theory. It is responsible for 
the modified hypcrscaling relation 2 — asp = {d—2)vBp. 
Based on the above considerations, it is reasonable 
to conjecture that, up to possible constants, A and g 
should be identified, as should x c — x and t uS. Thus, 

(r) 

apart from non-universal constants, G si ' is given by 
the one-point function (\&) in the supersymmetric the- 
ory (RW, which, by dimensional reduction, is the same as 
the one-point function (ip) in the Yang-Lee scaling the- 
ory (fllh. For d = 2 the gradient terms are absent, so 
G^l (x,g) — bi(gd/dx)lnZi where 



Zi = 



a (b2/fl)(-(* c -xW>+i^ 3 ), 



and b\ and 62 are non-universal constants. After rescal- 
ing the integration variable, this gives the main results 
(E|,p|) of RBG, together with the values 9 = \, (j) = | for 



the exponents. 

According to (g), the scaling function F(s) has sin- 
gularities at the zeros of the Airy function, which lie 
on the negative real axis. The closest to x = lies at 
x c (g) — x c + (2.388 . . .)(g/b 2 ) 2 ^ 3 , governing the asymp- 
totic behavior ~ x c (g)~ N of J2aPn A e ~ 9A as N — > 00, 
for fixed small g. This singularity is a simple pole, corre- 
sponding to the value asp = 3 — 9bp = 2. All this agrees 
with general crossover theory M that the scaling function 
should exhibit the critical singularities of the stable fixed 
point. 

We now discuss the generalization to area-weighted 
two-dimensional self-avoiding loops at higher order mul- 
ticritical points. The additional interactions between 
nearby portions of the self-avoiding loop will modify the 
parameters Up in (|7J) . As long as the truncation of terms 
leading to (g) remains valid, the dimensional reduction 
argument still applies with a modfied potential V, so that 
the rooted generating function is still given by the log- 
arithmic derivative of a generalized Airy integral of the 
form J c e~ v W / ' 9 dp . The obvious candidates for poten- 
tials which then yield multicritical behaviour in the limit 
g -> have the form |§ V(ip) = £* =1 Vjip j - ijj k+2 (the 
coefficient of ip k+1 is redundant, as it can be removed by 
a shift in ip.) Here v\ is linear in x c — x, and one can check 



that the other coefficients Vj should be positive deep in- 
side the single phase region. Repeating the above analy- 

(vj,g) = g(d/dvi)lnZ k 



sis then leads to the result G sing 

where 



Z k = 



-(EL^^'-^ +2 )/9 



df> 



with the contour C chosen to guarantee convergence. 
Comparing with the scaling form (ra) then gives the re- 
sults (0), together with 9 k = l/(fc + 2). In particular 
(R 2 ) ~ (N) 2u «, where v k = l/y x (k) = (k + 2)/2(fc + 1). 
It should be noted that, although these exponent values 
are based on extremizing a simple polynomial, they are 
not the same as those in Landau theory, in which the 
analogous potential would be X^=i v j4> 2 ^ + ^ 2 ^ fe+1 ^- 
Finally, the limit g — > may be taken in ( |l2[) , using the 

(r) 

saddle-point method, with the result that G si ' is sim- 
ply given by a zero of V'(ip). By considering the limit 
when all the Vj are large and positive, it may be shown 
that the correct zero in this single phase region is that on 
the real axis with the largest real part. The contour is 
to be run through this, locally parallel to the imaginary 
axis. Thus, near the k = 2 multicritical point, where 
V(ip) =v\i])-\- V2ip 2 — ip 4 i one finds the result (g), as the 

appropriate root of a cubic equation. This formula has 

(V) 
a number of interesting properties. At V2 = 0, G si ' be- 
haves like (x c — x) 1 / 3, . For fixed V2 > 0, the first singular- 
ity occurs not at s = 1 , but at s = — 1 , and this is a square 



'■-":<& 



(x c (v 2 )-x) 1/2 , where x c (v 2 )-x c (0) 



, 3/2 



all as expected on the basis of crossover theory. When 
V2 < 0, the corresponding scaling function is $ < (s) — 



.1/3 



((s 2 + l) 1 / 2 + s) l/6 - ((s 2 + l) 1 / 2 - s) 1/u , where now 
s = vi/(—V2) 3 ' 2 - The other two roots remain complex, 
and there is now a first-order transition at s « 1.6, when 
Re V has the same value at each root. Note that $ < , con- 
tinued into the phase coexistence region, is analytic near 
the real axis and does not exhibit a spinodal singularity. 
However, there is a puzzle, which our theory has in 
common with that of Saleur J7J, associated with the 
physical identification of this multicritical point. The 
generalisation of the Flory approximation flq] to the 
and higher-order multicritical points, at which the first k 
renormalised virial coefficients vanish, may be phrased as 
follows: consider a long loop of length N and linear size 
R. In the absence of interactions, its entropy may be esti- 
mated on the basis of a free random walk to be 0(R 2 /N). 
The mean density is 0(N/R 2 ), so that the interaction 
energy may be approximated by v k+ iN(N / ' R 2 ) k , where 
v k +\ is the first non-vanishing virial coefficient. Bal- 
ancing these two contributions to the free energy then 
gives R ~ N Vk with v k as above. This argument, com- 
bined with observation that at the fcth multicritical point 
there should be exactly k relevant parameters, suggests 
that k — 2 in our theory should be identified with the 
0— point, and that, because of the simple scaling associ- 



ated with the dimensionally reduced theory, the Flory re- 
sult is in fact exact. However, the value ^2 = 1, and the 
associated crossover exponent y 2 (2)/yi(2) = |, do not 
agree with the predictions of an exactly solvable model of 



Duplantier and Saleur |19|, for which the corresponding 



(12) values are j and |. While it might be argued that this 
model is somewhat special, extensive numerical studies 
of more generic lattice models appear to confirm these 
values J2(J . Even more strikingly, the supposedly correct 
value for v at the ©-point actually corresponds to k = 6 
in our theory (and that of Saleur) , and the crossover ex- 
ponent corresponds to a perturbation with j = 5. But 
this relevant perturbation should lead to a multicritical 
point with k = 3, not the usual one with k = 1 as ex- 
pected on the physical grounds for the ©-point. This 
same paradox was present in Saleur's theory [|7| and at 
present there seems to be no plausible resolution. It may 
well be that the coincidence of the ©-point exponents 
with k = 6 is merely that, and that the sequence of mul- 
ticritical points of self-avoiding loops implied by Saleur's 
and the present theory represent some completely differ- 
ent physics. Unfortunately, because of the truncations 
made in going from the original self-avoiding loop model 
to the generalised Airy integral, it is very difficult to say 
to what these other multicritical points might correspond 
physically. 

To summarize, the conjectured scaling function of 
Richard, Guttmann and Jensen M for area-weighted self- 
avoiding polygons has been shown to follow from physi- 
cal reasoning concerning the crossover to branched poly- 
mers, together with the dimensional reduction argument 
of Parisi and Sourlas [pi. Depending on ones point of 
view, the numerical confirmation of this formula found 
by RGJ could be taken as dramatic vindication of the 
dimensional reduction argument, beyond its simple pre- 
diction of the value of the entropic exponent asp = 2. 
The exact formula for the scaling function is in accor- 
dance with standard crossover theory |§|] , but it points to 
the importance of understanding all the singularities of 
the crossover scaling function, not just the physical ones, 
in building up the full scaling form. 

The formula proposed by RJG is just the first of a 
series of exact scaling functions describing higher-order 
multicritical points for self-avoiding loops weighted by 
their area. Moreover this approach enables one to recover 
exact results for scaling functions in the unweighted en- 
semble, and these have the form of algebraic functions. 
These are the first examples of exact but nontrivial scal- 
ing functions of more than one thermodynamic variable 
at isotropic critical points. 

The simple structure found here is analogous to that 
which appears in N — 2 supersymmetric theories in two 
dimensions, although in this case the supersymmetry is 
of a different nature. From that point of view the mul- 
ticritical points of self-avoiding loops correspond to the 
A k +\ series of simple singularities pl| of the potential V: 



it would be interesting to find analogs of the Dk+i series, 
and the exceptional cases. 
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